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Abstract 

Experimental data from a turbulent jet flow is analysed in terms of an additive, con¬ 
tinuous stochastic process where the usual time variable is replaced by the scale. We 
show that the energy transfer through scales is well described by a linear Langevin 
equation, and discuss the statistical properties of the corresponding random force in 
detail. We find that the autocorrelation function of the random force decays rapidly: 
the process is therefore Markov for scales larger than Kolmogorov’s dissipation scale 
r]. The corresponding autocorrelation scale is identified as the elementary step of 
the energy cascade. However, the probability distribution function of the random 
force is both non-Gaussian and weakly scale-dependent. 

PACS numbers: 47.27.Gs, 47.27.Jv, 05.40.-|-j 


1 Introduction 


Due to the molecular viscosity of the fluid, the kinetic energy of a flow must 
be dissipated. The degree of instability of a flow is usually quantihed by the 
Reynolds number Re, dehned as Re = UL/v where U, L and u are respectively 
the mean velocity, the largest scale of the flow, and the kinematic viscosity of 
the fluid. For large enough Reynolds number, the dissipation scale rj at which 
energy is dissipated becomes much smaller than the integral scale L at which 
it is injected: L ^ rj. The transport process through which energy injected 
at large scale is transported down to small scales is traditionally referred to 
as a “cascade”, following Richardson’s qualitative picture of turbulent eddies 
breaking down into smaller sub-eddies from the largest scales of motion down 
to the dissipative scales [1,2]. 
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Richardson’s ideas were first made quantitative by Kolmogorov [3]. Kolmogo¬ 
rov’s theory underlines the central role plaid by the mean energy dissipation 
rate (e). Its main assumptions are that, for sufficiently large Reynolds num¬ 
ber, (i) small scale turbulence (i.e. r L) is homogeneous and isotropic, 
and thus decoupled from the large-scale properties of the flow, which may be 
affected by, e.g., the specihc geometry of boundary conditions; {ii) for scales 
large compared to the dissipation scale (r 77), (e) is independent of v. 
Whenever ( 7 ) and {ii) hold, one expects small scale statistical properties of 
the cascade process to be universal, and determined only by (e) in the inertial 
range 77 r L. In addition, dimensional arguments lead to a simple scal¬ 
ing form for the probability distribution function of velocity increments. This 
theory becomes inconsistent if one allows the energy dissipation rate to fluctu¬ 
ate. However, experimental data shows unambiguously that the Kolmogorov 
picture is incomplete [1,2]: the energy dissipation rate e is a strongly fluc¬ 
tuating quantity. Its inhomogeneous, intermittent behaviour effectively links 
small and large-scale features of the flow. Further, the probability distribution 
function of velocity increments is not scale-invariant: its shape is Gaussian at 
large scale, yet develops long tails at smaller scale. 


The hrst phenomenological model to take into account the strong variability of 
the energy dissipation rate was proposed by Kolmogorov and Obukhov [4]. In 
this picture, emphasis is shifted from (e) to the local average integrated over 
a ball of radius r. The fluctuations of the random variable In e,. are assumed 
to be Gaussian: 


P(lner, r) 


1 ( (Ine,. - (lne,.))^\ 


( 1 ) 


Further, the variance A^(r) depends logarithmically on scale: 


J\^{r) = Aq + /iln(L/r), 


( 2 ) 


where p is a constant and Aq accounts for the large scale properties of the 
flow. Log-normal probability distribution functions (Eq. (1)) are known to be 
a good hrst-order approximation of the statistics of the energy transfer rate in 
fully-developed turbulent flows [5,6]. They arise naturally as solutions of mul¬ 
tiplicative stochastic processes for the variable [7,8], where interaction is 
generally assumed to be local in scale: the stochastic process must be Markov. 
However, this assumption is generally not verihed by experimental data [9]. 
Approaches related to the log-normal model (Eqs. (1-2)) have attracted re¬ 
newed attention recently [10,11], despite well-known difficulties linked to the 
scaling properties of and of the velocity increments 5vr- [8,12]. 

The statistical properties of isotropic and homogeneous turbulent flows have 
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been reconsidered recently within the context of continuous, additive, Markov 
stochastic processes [13,14], Analysis of experimental data from a turbulent 
flow has shown that the scale-dependence of probability distribution func¬ 
tions of velocity increments [13] and of the logarithm of the energy dissipation 
rate [14] may be described by Fokker-Planck equations, where the usual time 
variable is replaced by a monotonous function of scale r. The energy cas¬ 
cade corresponds to an Ornstein-Uhlenbeck process for the stochastic variable 
Ine^, since the drift and diffusion terms of the Fokker-Planck equation are 
respectively found linear and constant. This model is exactly solvable. The 
probability distribution function is found Gaussian: the statistics of are log¬ 
normal. However, the scale-dependence of the mean and variance differs from 
that considered in previously studied log-normal models [4,10]. The presence 
of a non-zero drift term implies that this new phenomenological approach is 
not equivalent to usual multiplicative models for the variable e^. 

In this work, we consider an alternative formulation of the model introduced 
in [14], Our analysis is based on the same velocity signal, recorded in a low- 
temperature turbulent jet flow at high Reynolds number Re = 20000. We 
aim at building a faithful, yet simple model of available experimental data, 
from a phenomenological perspective. Fokker-Planck equations, which govern 
the evolution of probability distribution functions, are equivalent mathemati¬ 
cally to stochastic differential equations with an additive noise term (Langevin 
equations, see [15]). The goal pursued here is to investigate the energy trans¬ 
fer process from the complementary viewpoint given by stochastic trajectories 
in scale. A detailed study of statistical quantihers which arise naturally in 
this new context allows to check the overall consistency of the approach. In 
particular, the validity of approximations made in [14] is investigated in detail. 

This article is organised as follows. The experimental set-up is discussed in 
Sec. 2, together with the specihcs of data acquisition and signal processing. 
Prescriptions used in the analysis of the experimental velocity signal are in¬ 
troduced in Sec. 3. The results previously obtained in [14] are summarised in 
Sec. 4. The stochastic cascade process is next examined in detail within the 
context of a Langevin description in Sec. 5. The Markov and Gaussian nature 
of the process are then discussed in Secs. 6 and 7 respectively. Sec. 8 is devoted 
to a summary and discussion of our results. 


2 A low-temperature jet flow 


The velocity signal analysed in this work was recorded in a low-temperature 
gaseous helium jet flow [16-18]. Experimental conditions (temperature T = 4.2 
K, pressure P = 1 bar) are set close to the liquid-gas transition: this ensures 
a very low kinematic viscosity, z/ ~ 2 x 10“® m^ s“^, and thus gives access to 
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high Reynolds numbers, up to Re ~ 10^ in this cell. 


After laminarisation, gaseous helium is injected through a contracting nozzle 
(diameter 2 mm) into a cylindrical cell (diameter 13 cm, height 30 cm) where 
the turbulent jet develops. A grid is set 17 cm downstream from the nozzle 
to avoid instability occurring on long time scales due to recirculation of the 
flow along the walls. Further downstream, gaseous helium is pumped out of 
the cryostat. Adjusting the microvalve inlet opening and the outlet pumping 
allows to vary not only the flow rate (mean velocity U), but also the operating 
pressure P in the cell. Since the product Pu is roughly constant, this original 
feature of the experiment allows to vary the Reynolds number over the range 
10^ < Re < 10^. Access to a large range of Reynolds number for a hxed 
geometry is the main beneht of this cryogenic experiment (see [18] for further 
details). 

The longitudinal component of the velocity is measured thanks to a cryogenic 
hot-wire anemometer [17] located 10 cm downstream from the nozzle on the 
axis of the jet, where the flow is fully turbulent. The effective spatial resolution 
is estimated to be 22pm. The wire is operated at constant resistance through 
a home-made 10 MHz lock-in ampliher. The signal is next digitalised using a 
12 bits A/D converter. We believe that the dominant source of experimental 
error is the resolution of the digitalization. 

What represents, from a physical point of view, the signal recorded by the hot 
wire? The velocity held is a function of space and time v{x, t). In a jet geom¬ 
etry, as in many experimental hows, the turbulent velocity held is advected 
relatively to the observer (anemometer). The “rate of velocity huctuations” 
is ususally dehned as the ratio of the standard deviation of the velocity over 
its mean value. A small value of this rate means that the how is ehectively 
“frozen”: Taylor’s frozen turbulence hypothesis is valid [1]. Rapid advection 
of the velocity held v{x,t) with respect to the (hxed) detector ensures that 
the signal records the spatial structure of the velocity held {v{x,t) ~ v{x), 
with X = {v)t), rather than its temporal evolution. In this jet, the rate of 
huctuation is 28%, in agreement with values typical of jet hows. The frozen 
turbulence hypothesis is now a source of error [19], in fact a bias on estimates 
of the spatial velocity held v{x). 

In the present work, this point is taken into account by resampling the signal as 
follows (see [20] for further details). The spatial coordinate a; of a record taken 
at time t is x{t) = /q v{t)dt. The bias on the velocity signal can be eliminated 
by considering instead the series v{x{t)). However, the interval between two 
successive points becomes irregularly distributed. One must also resample the 
signal regularly, in practice by linear interpolation. We checked that a more 
precise interpolation, involving the second order derivative, does not sensibly 
improve the regularity of sampling, and leaves our measurements unchanged. 
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Moreover, our conclusions do not seem to depend on the choice of a specihc 
correction method, since results remain unchanged when using, e.g., the local 
Taylor hypothesis advocated in [19]. 

For consistency, the velocity signal we analyse is the same as studied in [14], 
and was recorded at Re = 20000. The Reynolds number based on Taylor’s 
microscale is Rx = 341. The dissipation scale is estimated to rj = 20/rm 
according to Kolmogorov’s formula: 


V 



1/4 


( 3 ) 


with (e) defined as { (du/dx)^ ). The smallest structure of the flow is 
therefore resolved at this value of Re: this justifies our choice of this specihc 
sample. The integral scale, operationally dehned as the correlation length of 
the velocity signal, is evaluated to L ~ 1 cm (L ~ SOO?/). This numerical value 
indeed corresponds to the largest structure of the how, since the diameter of 
the jet is about 2.5 cm in the vicinity of the anemometer, where the measure¬ 
ment is done. Note however that the values of r] and L only have qualitative 
importance in our study. 


3 Prescriptions for data analysis 


As proposed by Obukhov [4], the dissipation averaged over a volume of size 
r may be used as a reasonable Ansatz for the local transfer rate of energy at 
scale r. Since available experimental data is often limited to the longitudinal 
velocity component v{x), it is customary to use the following one-dimensional 
surrogate of the dissipation er(x) at location x: 


er(x) 


x+rl2 . 2 

15 - r 


r 


x-r/2 


\dx') 


( 4 ) 


a dehnition which we endorse here. Eq. (4) is generally believed to capture 
most of the relevant physics (see [2] for up-to-date critical overviews). 

Even though the relevant physics is continuous, experimental data has been 
sampled before being recorded: the available data is thus discrete, with a 
sampling distance equal to d = 22/im (or d = t-tf])- In what follows, the unit 
of length is set to d unless otherwise noted. Values of the dissipation field 
er{x) are obtained thanks to a hrst-order hnite-difference approximation of 
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the derivative of the velocity held: 


15z/ 


x+int(r/2) —1 


er(X) = 


E 


x'=x—int{r/2) 


'v {x' + 5x') — v{x') 
5x' 


( 5 ) 


where 5x' is the (constant) discretisation step. The numerical value: 

5x' = 5, ( 6 ) 


is adopted in all calculations reported here. Note that with this choice, the 
huctuations of at r ~ 77 cannot be resolved since 5x' = 5.5r]. 

In the following, we will refer to both variables r and I as “scale”, where the 
new variable I is related to the physical scale r (as used in Eq. (4)) by : 

I = In (^) , (7) 

This dehnition is used for convenience since it leads to simpler analytic expres¬ 
sions. For the relevant scales Tj < r < L, the continuous variable I is positive 
and bounded from above: 

0 < / < In ^ ( 8 ) 

It is a monotonously decreasing function of the usual scale r. For clarity, scale 
will however be measured in units of the dissipation scale 77 in most graphs. 

The random process X{1) is dehned as: 

X(/) = lne„ (9) 

using Eqs. (4) and (7). Statistical averages, denoted by the symbol ( ), are 
evaluated over the whole data set (10^ points, or about 2 lO^L). For all sta¬ 
tistical quantihers presented here, this set is large enough to ensure statistical 
convergence, as conhrmed by comparisons with results obtained from smaller 
samples. 

For convenience, we will also consider the centered variable Y{1), dehned by: 

Y{l)=X{l)-{X{l)). (10) 


The probability distribution functions P{Y,r) of the variable Y{r) computed 
according to Eq. (5) are shown in Fig. 1. Despite rough agreement with the 
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parabolic shape expected for a Gaussian distribution in a lin-log plot, an 
obvious departure from the functional form (1) is the asymmetry of P{Y,r). 
This point has already been noted in previous experimental [5] and numerical 
work [6]. 



Y 


Fig. 1. Probability distribution functions P{Y, r) of the centered variable Y at scales 
r ranging from r = 25r/ to 400ry. The variance (Y(r)^) is a decreasing function of 
scale. Asymmetry of the histograms translates into a non-zero value of the skewness. 
For clarity, histograms were shifted by a factor 2 along the vertical axis. 

Evaluation of the scale-derivative of Y{1) will also be needed in the following. 
A hrst-order approximation will be used: 

^(0 = 4 (y(i + «)-K(/))i (11) 

where the step 51 is naturally related to the physical scale difference 5r by the 
relation: 


5l{r) = In 


' r + 5r' 


( 12 ) 


We use here the numerical value: 


5r = 4, 


(13) 


or 5r = AArj. The discretisation step is thus a decreasing function of physical 
scale r: 5l{r- = 25ri) ~ 0.17, 5l{r = lOOrj) ~ 0.04, 51{L) ~ 0.009. Our results 
are unchanged for other reasonable numerical values of the discretisation steps 
(6) and (13). This ensures that the continuous limit is well-controlled. 
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4 An Ornstein-Uhlenbeck process 


In this section, we wish to summarise the analysis made in [14] of the stochastic 
process X{1), as dehned above. 

A hrst important result is that the transition probability distribution function 
P {X',l'\X,l) = P (X(/')|X(/)) respects the Chapman-Kolmogorov equation 
to a very good approximation. Mathematically, the Chapman-Kolmogorov 
equation is a necessary (but not sufficient) condition for the process to be 
Markov. An equivalent (differential) formulation is the Kramers-Moyal ex¬ 
pansion [15], an evolution equation for the probability distribution function 
P(X, /) which reads: 

(9 oo an 

^P{X, 1) = y;(-i)”_ (D4X, 1)P(X, i)). (14) 


The Kramers-Moyal coefficients Dn{X,l) in Eq. (14) are dehned as the limit: 

P^X,/) = hmM„(X,/,/'), (15) 

where the functions M„ (X, /, /') are equal, for I ^ I' and up to a multiplicative 
factor, to the n-th order moment of the transition probability distribution 
function: 

M„ (X, /,/') = ^ i I (X' - X)- P(X', /'|X, /) dX'. (16) 

In [14], the coefficients Dn{X,l) are equated to numerical values found for 
M„(X,/,/' = I + SI), SI = 0.04, below which resolution becomes insufficient. 
According to this procedure, the second and fourth-order coefficients D2 and 
D 4 are found roughly constant, i.e. independent of both variables X and /, 
with P 4 < 0. 05 (^ 2 )^, and D 2 = 0.03 ±0.005. The inequality is interpreted by 
neglecting the fourth-order term in Eq. (14): P 4 (X, /) = 0. If P 4 = 0, Pawula’s 
theorem [15] then implies that the expansion (14) may be truncated at second 
order. In particular, D^{X,l) = 0, and Eq. (14) reduces to the Fokker-Planck 
equation: 

|P(V /) = -± {D,(X, l)P(X, l)) + -X (D,(X, l)P(X, /)), (17) 


where D 2 {X, /) = P is a diffusion constant. 

Using the same procedure, the hrst-order coefficient is found linear in X, with 



an approximately scale-independent slope 7 . Using at each scale the centered 
variable X — {X{1)), one may write: 

Di(X,/)=7(X-(X(/))) + F(/), (18) 

where F{1) is a (real) function of scale, and 7 = 0.21 ± 0.02. Eqs. (17)-(18) 
dehne an Ornstein-Uhlenbeck process with drift and diffusion coefficients 7 
and D [15]. The drift coefficient 7 , a real number, should not be confused 
with the drift term Di{X, /), a real function of variables X and /. Note that 7 
is found positive: the process does not relax to a stationary state for large /. 

In [14], Eq. (17) is integrated by assuming that energy does not fluctuate at 
the largest scale L, as expected in the ideal case of an ensemble of systems 
supplied with the same power. The corresponding initial condition reads: 

P(X,/= 0) = (5(X(0)-Inci) (19) 


The Eokker-Planck equation supplemented with Eq. (19) admits an exact, 
Gaussian solution: 


P(X,l) 


1 exp (_(£xk£M)!i 
A{l)V^ 2A2(0 j ■ 


( 20 ) 


Similarly to Kolmogorov and Obukhov’s model [4], the dissipation exhibits 
log-normal fluctuations. 

Since energy is conserved through the cascade (i.e. Vr G [ 7 , L], (e^) = (e)), the 
mean and variance of the Gaussian probability distribution function (20) are 
simply related by: 

{X{l)) = {Xm-A\l)/2, (21) 


in good agreement with experimental data. According to Eq. (19), we postulate 
that (X(0)) = X(0) = Ine^. The variance A^(/) of the Gaussian solution (20) 
is a monotonously decreasing function of scale r: 



A consistency condition is: 

F(l) = 4(A'(/)) = (23) 
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in good agreement with direct measurements of the first-order Kramers-Moyal 
coefficient -Di(X, /). Note that Eq. (2) is recovered when the drift coefficient 7 
cancels, with /r = 2D and Aq = 0. For small enough scales {r L), Eq. (22) 
reduces to a power-law dependence, in agreement with Castaing’s model. In 
that context, 7 was observed to be inversely proportionnal to log(i?e) [10,16]. 

One may summarise the main result of [14] as follows: an Ornstein-Uhlenbeck 
process with positive drift coefficient (Eqs. (17)-(18)) gives a simple picture 
of the scale dependence of energy transfer statistics of a turbulent flow. The 
validity of this description relies on two assumptions: first that the process 
is Markov, second that the Kramers-Moyal coefficients of order strictly larger 
than 2 may be neglected, in particular that D^{X,l) and D/^{X,l) may be 
set to zero. The analysis made in [14] suggests that both assumptions are 
reasonably well-supported by experimental data. However, one may already 
note that the non-zero, negative skewness of P(X, 1) observed in Fig. 1 is not 
taken into account by the Fokker-Planck model. In what follows, we adopt 
the complementary perspective provided by a formally equivalent Langevin 
description in order to better assess the validity of these approximations. 


5 A Langevin equation 


We will study in this section the energy cascade process from the perspec¬ 
tive given by scale-dependent random trajectories of the centered stochastic 
variable Y{1). A Langevin equation formally equivalent to the Fokker-Planck 
equation (17) is first introduced in Sec. 5.1. We proceed to measure directly 
the drift and diffusion coefficients 7 and D (Secs. 5.2) by methods which do not 
involve the computation of Kramers-Moyal coefficients. Predictions involving 
the two-point correlators of Y{1) are next confronted to experimental data in 
Sec. 5.3. A definition of the random force is finally introduced in Sec. 5.4, 
thanks to proper discretisation of the Langevin equation. 


5.1 The equivalent Langevin description 

Eqs. (17-19) and (23) are mathematically equivalent to the stochastic differ¬ 
ential equation [15]: 

ny 

—d) = jY{l) + V2B «pp(0 (24) 

with initial condition K(0) = 0 (note that F{1) = d{X{l))/dl). When strict 
equivalence with the Fokker-Planck equation (17) holds, the random force 
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.^Fp(0 is a stationary, Gaussian, white noise and the stochastic process is 
Markov. The two-point correlation function therefore reads: 

('Cfp(0 = S{1 - /'), (25) 


and the probability distribution function is a scale-independent Gaussian: 

with mean (,^fp(0 ) = 0 variance (.Cfp(O^) = 1 - 

The linear equation (24) is exactly solvable. Its solution reads: 

i 

Y{1) = F(0) + J ^{l') dl'. (27) 

0 

Two-point statistics of the process Y{1) can also be calculated when the ran¬ 
dom force is delta-correlated (Eq. (25)). Assuming that the initial distribution 
P(E, 0) admits a hnite width (y(0)^), one obtains: 

(F(/)2) = (y(o)2) + £ . (28) 


Expression (22) is recovered when {Y (0)^) = 0, as expected. Further, we dehne 
the normalised scale autocorrelation function of the stochastic process Y{1) by 
the relation: 


Cy{1,AI) 


{Y{1) Y{l + Al)) 

{ym 


(29) 


where Al is a positive scale difference. A straightforward calculation leads to: 


lnGy(/,A/) =7 A/. 


(30) 


Note that expressions (28) and (30) are obtained independently of the func¬ 
tional form of the probability distribution function of the random force. 


5.2 Drift and diffusion coefficients 


Assuming that the energy cascade process can be described by a Langevin 
equation such as (24), one would hrst of all like to check whether this equation 
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<AY/At I Y > 


is indeed linear. Since (^(/)) = 0, a straightforward consequence of Eq. (24) is: 


(y(/)|r(/))=7r(/), (31) 


where Y (/) denotes the derivative of Y (/) with respect to scale. It is easy to 
see that nonlinear functions of Y (/) in the right hand side of Eq. (24) would 
also contribute to the the conditional average {Y{l)\Y{l)). Direct measures 
of {Y{l)\Y{l)) are presented in Fig. 2, where the prescriptions detailed in 
Sec. 3 were used. Within statistical error, the conditional average {Y{l)\Y{l)) 
is indeed proportional to Y (/), even though a small, possibly cubic contribution 
appears at small scales for (rare) large negative values of the variable Y (/). This 
observation conhrms the linearity of the relevant Langevin equation. However, 
the corresponding slope 7 seems to be an increasing function of scale. 

(a) (b) 




Fig. 2. Measure of the drift coefficient 7 (/). Graph (a): plot of the conditional average 
{Y{l)\Y{l)) vs. Y{1), estimated at scales r = 257,1007 and 4 OO 7 . Straight lines of 
slope 0.20, 0.24 and 0.30 are drawn to guide the eye. Graph (b): plot of 7 (r), as 
defined in Eq. (33), versus scale r/ 7 . Statistical error is of the order of a few percent 
in relative value. Inset; log-lin plot of the same data (I = ln(L/r)). The straight line 
corresponds to 'y{l) = 0.32 — 0.05 1. 

Using Eq. (31) to obtain quantitatively accurate estimates of 7 (/) is rather 
costly numerically. We will turn to a simpler method, which turns out to yield 
consistent estimates. Implicit in Eq. (24) is the assumption that the random 
force ^Fp(0 is statistically independent from the stochastic variable Y{1) at 
all scales /: 


{Y{1) eFp(/)) = 0. 


(32) 
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Upon multiplying both sides of the Langevin equation (24) by Y (/), ensemble¬ 
averaging, and using Eq. (32), one easily obtains: 


7(0 


(H) yji)) 

{ym 


(33) 


where the drift coefficient 7 (/) is a priori a function of scale. Fig. 2.b shows 
that 7 (/), as estimated thanks to Eq. (33), is indeed a slowly varying function 
of scale. The order of magnitude is however the same as that of the (constant) 
value advocated in [14] (0.21±0.02). This method appears to be more sensitive 
than direct calculations of the Kramers-Moyal coefficients. 


We would like to emphasise that: (i) Eqs. (31) and (33) yield mutually con¬ 
sistent estimates of 7 (/); {ii) the normalised cross-correlation function Uyg(r), 
dehned as: 


Crdr) 


{Y(r) «r)) 


(34) 


is indeed close to zero (in practice of the order of 10 “^) when the random force 
^{l) is calculated thanks to Eq. (43) for constant drift and diffusion coefficients, 
e.g. 7 (/) = 7 = 0.21 and D{1) = D = 0.03. 

As shown in the inset of Fig. 2 , the scale dependence of the drift coefficient is 
well described by a linear function of I over the whole range of relevant scales. 
Introducing a factor 2 for convenience, we write: 


7(0 = 70 - 27i(. 


(35) 


Our estimates of the constants 70 and 71 are 70 = 0.32 ± 0.03 and 71 = 
0.025 ± 0.003. The error bars take into account both statistical error and 
the uncertainty deriving from the existence of other possible choices of the 
discretisation steps Sr and Sx'. Note that the drift coefficient is positive for all 
physically relevant scales: Eq. (35) yields 7 (r = 7 ) = 0.01 > 0. 

Assume now that the second-order Kramers-Moyal coefficient D 2 {X, 1) is inde¬ 
pendent of X, and reduces to a (possibly scale-dependent) diffusion coefficient 
D{1). According to Eq. (24), this coefficient may then be dehned as: 

2fl(0 = ((^(0-7(/)r(/)) ). (36) 

where the variance of the random force is set to unity: = 1. Using 

Eq. (36), where 7 (/) is evaluated thanks to Eq. (33), we hnd a hnite, scale- 
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dependent, monotonously decreasing coefficient diffusion D{r) (as a function 
of the physical scale r, see Fig. 3). The order of magnitude is consistent with 
that of the constant value advocated in [14] [D = 0.03 ± 0.01). An excellent 
£t of the data is given by the expression: 


D{1) = D, e^ 


51 


(37) 


with Do = 0.01 ± 0.005 and S = 0.40 ± 0.01 (see the inset of Fig. 3), except at 
scales r close to the dissipation scale rj. 



rA] 


Fig. 3. Graph of the diffusion coefficient D(l) (as defined in Eq. (36)) versus 
scale r/rj. Statistical errors are smaller than the width of the symbols used. In¬ 
set: same data in log-log scale. The straight line corresponds to the scaling law 
D(r) = 0.01(L/rf-^. 


We have shown in this section that the drift and diffusion coefficients, when 
measured according to Eqs. (33) and (36), are in fact not independent of the 
scale /. Since this scale-dependence cannot be attributed to, say, statistical 
errors, it becomes essential to understand to what extent the coefficients j(l) 
and D(l) may be approximated by constant coefficients, as has been advocated 
in [14]. This is the goal of the next section, where predictions concerning the 
two-point statistics of y(l) will be confronted to experimental data for both 
scale-dependent and constant coefficients. 


14 




5.3 Two predictions 


Instead of Eq. (24), we now use the following Langevin equation as our starting 
point: 


— (/)=7(/)n0 + y^e(0, (38) 

where the functions 7 (/) and D(l) are given by Eqs. (35) and (37) respectively. 
The solution of Eq. (38) reads: 


Y{1) = 




f _ I- ' 

y(0) + f 


(39) 


Assuming that the random force ^{l) is delta-correlated, we obtain the follow¬ 
ing expression of the variance of Y (/): 


(K(/)^) = I {y(0)2) + 2D„J e2(«-To)r+27ii"j(' 


(40) 


for an initial width (y(0)^). The normalised scale autocorrelation function of 
the stochastic process Y{1) then reads: 

In C'y(/, A/) = ( 7 o- 27 i/)A/- 7 i(A/) 2 . (41) 


In Fig. 4.a, we compare the scale dependence of (Y{lY), as obtained from 
experimental data, with predictions (28) and (40). The initial value is set to 
(F(0)^) = 0.23. We use the following numerical values: 7 = 0.21 (Eq. (28)), 
7 o = 0.32, 7 i = 0.025, Dq = 0.01, 5 = 0.40 (Eq. (40)). We hnd that Eq. (40), 
not Eq. (28), hts experimental data extremely well. A small deviation is ob¬ 
served only for scales smaller than r = 25^ (/ > 3.0), where the diffusion co¬ 
efficient D{r) deviates from the scaling law (37) (see Fig. 3). This first result 
suggests that approximating the drift and diffusion coefficients to a constant 
value is inappropriate. 

Fig. 4.b shows the measured scale autocorrelation functions Cy(/, A/). Again, 
the linear growth predicted by Eq. (30) is not observed. However, the following 
features of the scale-dependence of Cy(/, A/) are reproduced by Eq. (41) at a 
qualitative level. The autocorrelation function is a function of both I and A/. 
For all (fixed) values of /, its logarithm first grows with A/ before reaching a 
maximum and eventually decaying. The tangent at A/ = 0 is smaller than 70 = 
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0.32. However, quantitative agreement is still lacking, for instance concerning 
the location of the maximum of C'y(/,A/), as well as the exact dependence 
in A/. We will show in Sec. 6 that, as one would now expect, a ^-function 
is indeed a good approximation of the autocorrelation function of the noise 
C^{l,Al), as shown by comparing the correlation scale of to the evolution 
scale of Cy- We therefore believe that the quantitative discrepancy observed 
between the prediction (41) and experimental data should be attributed to 
the approximation made when htting 7 (r) by a logarithmic function of scale. 
In particular, the oscillations apparent in the inset of Fig. 2 for large I do not 
seem to be statistical fluctuations, and should therefore be taken into account 
when computing Cy(/, A/). 

(a) (b) 




Fig. 4. Graph (a): variance of Y(l). The statistical error on experimental data 
is smaller than the width of symbols. The prediction obtained in the case of 
scale-dependent coefficients (Eq. (40), dashed line) is indistinguishable from ex¬ 
perimental data in the interval 0 < I < 3. See text for the numerical values used 
in the calculation of {¥{1)'^). Graph (b): lin-log plot of the scale autocorrelation 
functions Cy{1,AI) vs. the scale difference Ah for scales I ranging from I = 25r/ 
to I = 400^. A straight line with slope 70 is also shown for comparison with the 
behaviour expected in the vicinity of AZ = 0 in the case of scale-dependent drift 
and diffusion coefficients (Eq. (41)). 


In conclusion, 7 (/) and D{1) should not be approximated to a constant value. 
Further, the functional dependence given by Eqs. (35) and (37) is consistent 
with experimentally obtained two-point correlators of the process Y{1), quali¬ 
tatively only for the autocorrelation function Cy{1, A/), and quantitatively for 
the variance (Y{lY). 
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5-4 Discretisation scheme 


In the following, we will consider the Langevin equation (38), where 7 (/) and 
D{1) are defined by Eqs. (33) and (36) respectively. Provided that the noise 
term S,{t) is Gaussian, stationary, and delta-correlated, this equation is for¬ 
mally equivalent to a Fokker-Planck equation such as Eq. (17), with scale- 
dependent drift and diffusion coefficients. As before, this Fokker-Planck equa¬ 
tion admits a Gaussian solution, the variance of which is given by Eq. (40. 

The stochastic differential equation (38) may be discretised according to the 
following hrst-order scheme [15]: 

Y{1 + 61) - Y{1) = 7(/) Y{1) 61 + ^2 D{1) 61 ^(/). (42) 


In the following, we will use Eq. (42) as an operational dehnition of the driving 
random force: 

(r + ■*') - r (0 - 7(() Y(l) 61). (43) 

By construction, ^{l) is a centered stochastic variable at all scales /: {01)) = 0, 
with unit variance (■^(0^) = 1- 

In Secs. 6 and 7, we will respectively test the Markov and Gaussian nature of 
the stochastic process defined by Eq. (38). 


6 A Markov process 


In this section, we further investigate the validity of Eq. (38) by turning to the 
dynamics of the process, and in particular to two-point scale autocorrelation 
functions. 

Due to the viscosity of the fluid, the velocity held of a turbulent how remains 
diherentiable. The energy cascade cannot be perfectly represented by a Markov 
process: the process 0 ^) must also be diherentiable, and can therefore not be 6 - 
correlated. Its (normalised) two-point scale autocorrelation function is dehned 
as: 


cohi') = {01)01'))- 


(44) 


Fig. 5.a shows that the scale dependence of G^(/, I') cannot be reduced to the 
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CJl, Al) 


form C^{1' — 1). The process is thus instationary in terms of the scale variable 
1. One may however write without loss of generality: 

Q(/,/0=Q(/,A/), (45) 


where A/ = /' — / is set positive by convention. Further, correlations of the 
random force differ from the ideal case (Eq. (25)): the autocorrelation scale 
r(/) is non-zero, its numerical value depends on scale 1. Even though the scale 
autocorrelation function decays approximately exponentially over roughly one 
decade, we choose to evaluate the characteristic correlation scale r(/) from the 
expression : 


t{1) = J C^{l,Al)dAl. (46) 

0 


We checked that this particular choice does not affect our conclusions. 

The autocorrelation scale r(/) is a quantitative measure of the departure from 
the Markov approximation. At a given scale /, the scale r(/) marks the limit 
below which the stochastic cascade process becomes smooth. In this sense, the 
(measurable) scale r(/) is the elementary step of the energy cascade. 

(a) (b) 




Fig. 5. Graph (a): scale autocorrelation function C^{l,Al) of the random force ^(Z) 
for scales r = Lexp(—^) ranging between r = 25r] and r = 400r]. The statistical 
error is smaller than the width of symbols. Inset: lin-log plot of the same curves. 
Graph(b): scale-dependence of the autocorrelation scale r(r). Inset: same graph as 
(a), represented as a function of the physical scale r. 

Eq. (46) yields numerical values ranging between r = 0.37 when r = 25r] and 
r = 0.03 when r = L = 500?7. Fig. 5.b shows that the scale-dependence of r(/) 
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may be fitted by the expression: 


r(r) = In 




(47) 


for a numerical value of tq close to 11 rj. In other words, the instationarity 
observed at the level of autocorrelation functions (Fig. 5.a) is mostly due to 
the change of variable I = ln{L/r). This is conhrmed by the insert of Fig. 5.b, 
where the autocorrelation functions are plotted with respect to the physical 
scale difference Ar: the first decade of decay is characterised at all scales by 
a slope roughly equal to 1/ro. In terms of the scale variable r, the elementary 
cascade step may simply be dehned as ro, of the order of Kolmogorov’s scale 
V- 

The characteristic evolution scale of the random process Y (/) at scale I is given 
by the inverse drift coefficient l/7(/) (cf. Eq. (38)). The product 7 (r)r(r) is 
a monotonically decreasing function of scale r, with 7(25?7)r(25?7) = 0.06 and 
7 (L)r(L) = 0.01. In the range of scales 25?7 < r < 500?7 at least, one hnds 
that: 

r(l) « -I; 



by more than one order of magnitude. Although the correlation scale r(/) is 
non-zero, the random process Y{1) is therefore effectively Markov. Moreover, 
the product 7 (r)r(r), as dehned by Eqs. (35) and (47), admits an absolute 
maximum close to 0.1 in the vicinity of r ~ 5r]. Assuming that Eqs. (35) 
and (47) faithfully describe the scale dependence of 7 (/) and r(/) down to 
the smallest scales, this implies that the energy cascade is always a Markov 
stochastic process for scales larger than rj. 


7 A non-gaussian random force 


In this section, we investigate the validity of Eq. (26) by evaluating the prob¬ 
ability distribution function P{^, 1) of the random force ^ at scale /. 


7.1 A typical realisation of the process f{l) 


A typical realisation of the random processes Y(1) and f{l) is presented in 
Fig. 6. The dynamics of Y (/) is generally dominated by the deterministic part 
of Eq. (38). Long periods of quasi-exponential growth controlled by the drift 
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coefficient 7 > 0 are perturbed by small deviations due to the random force 
^{l). This regime is interrupted infrequently by sharp drops, corresponding to 
large negative excursions of the random force. This suggests that the probabil¬ 
ity distribution function of ^ is strongly skewed, in contrast with the Gaussian 
prediction. Accessible values of also seem to be bounded from above. 



15 I-^^^^^^^-1 

0.0 1.0 2.0 3.0 4.0 


I = ln(Ur) 


Fig. 6. A typical random trajectory of the stochastic processes Y{1) and ^{l). Note 
that Y (0) is not equal to zero. 


1.2 Probability distribution function of the random force 


As shown in Fig. 7, the probability distribution function P{f, r) of the random 
force f at scale r is markedly different from the Gaussian prediction (Eq. (26)): 
the skewness factor of P{f, r) is indeed negative. The shape of P(^, r) depends 
on scale r: this is consistent with the instationary behaviour of autocorrela¬ 
tion functions (cf. Sec. 6 ). Even though Gaussian-like, rapid decay of P{f,r) 
is observed for .^ > 0 , the probability of large negative deviations is much 
larger than that predicted for a Gaussian process. This long tail can be htted 
with reasonable accuracy either by a log-normal or by a stretched-exponential 
functional form. 

The Gaussian-like decay of P{^) for positive values of the random force may be 
understood as follows. The dehnition of as a sum of non-negative quantities 
(Eq. (4)) implies that: 


r' Cr' < r Cr, 'i r' < r. 


(49) 
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The process X{1) can only adopt values such that: 




(50) 


Taking the appropriate limit, one obtains: 

y(/)<i-d(x(/)). (51) 

The random force ^ is dehned as the difference of the two random variables Y 
and Y, where Y admits an upper bound and the fluctuations of Y are nearly 
Gaussian: the tail of 1) is thus expected to be Gaussian-like for positive 
values of the random force 



Fig. 7. Probability distribution functions P{^,r) of the random force ^ for scales r 
ranging from r = 25r] to r = 400?] in lin-lin (left-hand side) and lin-log (right-hand 
side) plots. The dashed line is drawn for comparison with the Gaussian prediction 
derived from the Fokker-Planck model. 

At a qualitative level, the asymmetry of the probability distribution function 
of ^ is readily understood if one remembers that P{Y, r) is also asymmetrical, 
as seen in Fig. 1. The solution of a linear stochastic evolution equation with 
Gaussian random driving is symmetrical. Gonversely, the skewness of the vari¬ 
able Y (/) points at the necessity of corrections to the Fokker-Planck model. 
The dominant one is expected to stem from the third-order coefficient D^. 
Non-zero n-th order coefficients M„ translate into deviations of the n-point 
correlators of the random force ^(/) from their Gaussian form. In particular, 
a non-zero value of D^, translates into a non-zero skewness, as observed here. 
The asymmetry of P(F, r) is weak, as quantified by the small value of the 
parameter |Zi) 3 (/)|/Zi)(/)^/^, which we found smaller than 10“^ at all scales. 
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Intermittency of fully-developed turbulent flows is often characterised by the 
presence of large excursions of the (positive) local dissipation rate e{x) = 
lbv{dv/dxY far above its mean value. We believe that these rare events also 
correspond to large negative deviations of the process ^{l). The dissipation rate 
er{x) at location x is dehned as the spatial average of e{x) over an interval of 
length r (Eq. (4)). The location x being hxed, er{x) is a monotonously decreas¬ 
ing function of scale r when the integral of e{x) over an interval of length r 
grows more slowly than r, i.e. when fluctuations of e{x) are weak enough. This 
case corresponds to the periods of quasi-exponential growth of Y (/) in Fig.6. 
Assume now that a sharp increase of e{x) takes place close to x* = x + r*/2. 
Close to r*, the value of er{x) will increase sharply: er*+&r{,x) 3> er*-&r{,x), cor¬ 
responding to a large positive value of the scale derivative der{x) /dr\r=r* ■ Since 
the logarithmic scale I is a monotonously decreasing function of the physical 
scale r, this shows that intermittent bursts of e{x) do indeed correspond to 
sharp and localised drops of Y{1). In other words, the asymmetry of P{^,r) 
expresses the intermittent character of the local dissipation held. 

However, the dehnition of as the averaged dissipation rate over a segment of 
length r is to some extent arbitrary. The amplitude of intermittent events, as 
quantihed by the value of ^{l*) close to a large deviation of the local dissipa¬ 
tion rate e{x*), may depend sharply on the precise dehnition of the averaging 
procedure which leads to er{x). In this sense, we conjecture that the precise 
numerical value of, e.g., the skewness coefficient of P{^,r) is in fact not rele¬ 
vant to the description of the physical cascade process, since it may depend on 
arbitrary features of the analysis method, such as the choice of an averaging 
window in Eq. (4). 


8 Discussion 


Following [14], we have found additional evidence that the energy cascade of 
a turbulent how at high Reynolds number is well described by a continuous 
stochastic process of Ornstein-Uhlenbeck type. We also clarihed the validity 
of a number of approximations made in [14]. 

Dehning the stochastic process Y{1) as the centered variable Ine^ — (Ine^) at 
scale / = ln(L/r), the relevant Langevin equation is linear, and reads: 

dY ! - 

aO- (52) 

We found that the scale-dependence of the drift and dihusion coefficients 7 (/) 
and D(/), as evaluated directly from experimental data by methods indepen¬ 
dent from those used in [14], is well described by the following functional 
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forms: 


7(0 = 7o - 7i 0 
D{l)=Do 


(53) 


where 70 = 0.32 ± 0.03, 71 = 0.025 ± 0.003, Dq = 0.01 ± 0.005, and 5 = 
0.40 ± 0.01 are positive constants. The (exact) solution of Eqs. (52)-(53) is in 
good agreement with experimental data at the level of two-point correlators 
of the stochastic variable Y{1). Further, we have shown that the drift and 
diffusion coefficients cannot be approximated by constant values if one wants 
to preserve this agreement. 

The main characteristics of the random force ^{l) have been determined di¬ 
rectly from experimental data. We have checked that the process is Markov 
at all scales /, since scale autocorrelation functions of the random force de¬ 
cay rapidly on a characteristic autocorrelation scale r(/) much smaller than 
the typical evolution scale l/ 7 (/) of the process. This is perhaps our most 
important result: the Langevin equation (52) dehnes a driving random force 
^{l) which respects the Markov hypothesis. Previous work [14] only provided 
evidence for the validity of a necessary condition for the additive process (38) 
to be Markov. The autocorrelation scale r(/) is the elementary step of the 
cascade process. It depends on the physical scale r as r(/) = ln((r -|-ro)/r), 
where tq is of the order of Kolmogorov’s scale 7 . 

Finally, our analysis of experimental data shows that the probability distribu¬ 
tion function of the random force ^{l) is strongly non-Gaussian, and exhibits 
a long tail for negative noise (see Fig. 1). These rare, but intense deviations 
are correlated with large positive values of the scale derivative der/dr. For a 
linear equation such as (38), the non-zero skewness of P{^,1) translates into 
an asymmetrical solution P{Y,l), in qualitative agreement with observation. 
This deviation from the simple Ornstein-Uhlenbeck picture is equivalent to a 
correction to log-normal statistics for the energy dissipation rate e^. Within 
the Fokker-Planck description, the corresponding leading-order correction is 
the third-order term of the Kramers-Moyal expansion (Fq. (14)), with a non¬ 
zero, scale-dependent coefficient This correction is weak, as quantihed 

by the value of the ratio \D 3 {l)\/D^P, smaller than 10“^ at all scales. This 
correction is made at the expense of the solvability of our model. 

Two comments are in order. First, this stochastic description of the energy 
cascade (including the Markov property) applies to all scales, from the integral 
scale L down to the dissipation scale 7 . The distinction commonly observed 
between dissipative and inertial sub-range statistics of, e.g., turbulent velocity 
and passive scalar density helds [ 1 ], appears to be irrelevant in the case of 
the energy dissipation rate. Second, all known multiplicative models of the 
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cascade use, at least to our knowledge, a scale-independent discretisation step 
such as In 2. Although natural from a mathematical viewpoint, this choice 
overlooks the scale-dependence of r(/) observed here, and may thus not be 
appropriate on physical grounds. The choice of the discretization step of a 
discrete stochastic model of the energy cascade is not meaningless. 

We would now like to emphasise that observations reported in this work stem 
from the analysis of one turbulent velocity signal, recorded in a particular 
realisation of a jet flow at a given value of the Reynolds number. One would 
of course like to know to what extent the Langevin equation (52) is a universal 
description of the energy cascade. Preliminary studies have already shown that 
the statistical properties of the processes ^{l) and Y (/) are qualitatively similar 
to those discussed here for other values of the Reynolds number [21]. Further, a 
number of adimensional quantihers of the energy cascade have been introduced 
at a rather formal level (7(/), D{1), etc.). One would like to understand their 
physical meaning, as indicated by, e.g., a possible dependence on the physical 
properties of the fluid and the Reynolds number, and in particular to elucidate 
their behaviour in the limit of large Re. Conversely, it is also essential to check 
whether some of our results, such as the numerical value of, say, the skewness of 
r), may not be artefacts of the data processing method. These important 
points are currently being considered, thanks to the analysis of the velocity 
helds of other turbulent flows. 

Interestingly, the existence of a stochastic evolution equation in scale allows, at 
least in principle, to reconstruct by simple integration the small scale statistics 
of the dissipation held from spatial huctuations observed at larger scale. We 
therefore believe that this description of the energy cascade process may prove 
useful, as a model of the small scale huctuations, in the context of large eddy 
simulations where only the larger scales of a how are usually resolved. 
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